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Abstract. Our efficient numerical approach has been appHed to model- 
ing the asymmetric circumstellar dust disk around (5 Pictoris as observed 
with the HST/STIS. We present a new model on the origin of the warp- 
ing of the P Pic disk. We suggest that the observed warp is formed by 
the gravitational influence of a planet with a mass of about 10 masses 
of Earth, at a distance of 70 AU, and a small inclination (~ 2.5°) of 
the planetary orbit to the main dust disk. Results of our modeling are 
compared with the STIS observations. 



Introduction 



In circumstellar disks, the major sources of dust are thought to be minor bodies, 
like comets. Both the stellar radiation drag (Poynting- Robertson drag) and 
stellar wind drag tend to induce dust inflow toward the star. As the dust passes 
by planets in its infall, it interacts with them, for example, by accumulating in 
planetary resonances (Liou &: Zook 1999, Ozernoy et al. 2000). The inner region 
of the P Pictoris dusty disk is tilted by a few degrees with respect to the outer 
disk (Burrows et al. 1996, Mouillet et al. 1997, Heap et al. 2000 (for references, 
see the last paper). Heap et al. (2000) concluded that the size and shape of the 
warp favors the presence of a planet (s) in a slightly inclined orbit. 



1 



APS Conf. Ser. Style 



2 



2. Modelling /? Pictoris Disk with Embedded Planet 

The dynamics of interplanetary particles are determined by several effects, in- 
cluding stellar radiation pressure; the Poynting-Robertson (P-R) and stellar 
wind drags; resonance effects by planets; and gravitational encounters with plan- 
ets. 

To simulate the gravitational dynamics of comets and nonconservative dy- 
namics of dust particles around a star, we use an implicit second-order integrator 
(Taidakova & Gorkavyi, 1999). In the absence of a planet, the surface density of 
dust between the star and the pericenter of dust sources is constant, while the 
number density of dust varies as r~^, r being heliocentric distance (Gorkavyi et 
al. 1997). The presence of a planet induces dramatic changes in the cometary 
and dust distribution. This is illustrated in Fig. 1, which shows a representative 
trajectory of dust particles. 
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SEMIMAJOR AXIS, AU 

Figure 1. Representative trajectory of a dust particle (/3 = 0.1) 
in the disk of /3 Pic. P-R drag and a slow precession of the particle's 
orbit around the orbital plane of the planet (a=70 AU; i = 2.5°) are 
the dominating dynamical processes during its initial inflow (220-83 
AU). The particle's inclination oscillates from to 2ipianet = 5°. The 
vertical line marks the resonant capture into the 4:3 resonance. The 
next 'jump' to ~ 17° is due to scattering of the particle by the planet. 
The final part of trajectory (< 60AU) marks the particle's drift. After 
strong interactions with the planet, dust particles (and comets) moves 
symmetrically about the new plane of the planetary orbit. 
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STIS OBSERVATIONS 
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Figure 2. Warp of the inner part of the /3 Pictoris disk. Upper panel: 
STIS/CCD coronagraphic images of the /3 Pictoris disk normahzed to 
the maximum flux with the vertical scale expanded by 4X (Heap et al. 
2000). Bottom panel: model disk with Mpi = 1OM0, orbital radius 
r=70 AU, and inclination i = 2.5°. The lifetime of dust particles due to 
collisions could be as small as 100 years. The stellar radiation pressure 
is negligible (dust particles are large enough). This warp is stationary, 
in contrast to the model of Larwood and Papaloizou (1997). 

Two algorithms for modeling the f3 Pictoris disk were considered. 
Algorithm I: 

1. Input an initial symmetric cometary disk. 

2. Input the mass and orbital radius of the planet on an orbit inclined to 
the initial cometary disk. 

3. Simulate the asymmetric dust distribution on a time scale of 1.4 x 10® 
yrs. Inclination of particles orbits to the plane of the initial cometary disk and 
asymmetry of dust is the result of precession of drifting particles around the 
planet's orbital plane (see Fig. 1). 

4. Simulate the dust-scattered light distribution. 
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Algorithm II: 

1. Input an initial symmetric cometary disk. 

2. Input the mass and orbital radius of the planet on an orbit inclined to 
the initial cometary disk. 

3. Determine the asymmetric distribution of comets due to orbital preces- 
sion and gravitational scattering by the planet (these processes can be realized 
without any drifts) after 2.8 x 10^ yrs. 

4. Simulation of the asymmetric dust distribution for short life times of 
particles. Asymmetry of dust is mainly a result of an asymmetry in the initial 
cometary distribution. 

5. Simulate the dust-scattered light distribution (see Fig. 2). 

Both alternatives produce a similar warp, but the second alternative seems to 
be more preferable due to: 

• a more realistic (shorter) lifetime of particles; 

• accounting for gravitational interactions between comets and the planet. 

The STIS images show the disk to be warped at small distances to the star in a 

sense that the inner disk is tilted by 4.6° with respect to the outer disk (Heap 
et al. 2000). See Fig. 2 to compare the STIS observations with our modeling. 

3. Conclusions 

Asymmetric structures in the circumstellar disk indicate at least one planet 
embedded in the disk of /3 Pictoris. 
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